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ABSTRACT 

Optical spectroscopic observations of the companion star (type G 8 IV to K4III) in the microquasar system XTE 
J1550-564 reveal a radial velocity curve with a best fitting spectroscopic period of Psp = 1.552 ± 0.010 days and a 
semiamplitude of JQ = 349 ± 12 km s“'. The optical mass function is/(M) = 6.86 ±0.71 Mq (Itr). We tentatively 
measure the rotational velocity of the companion star to be Uot sin; = 90 ± 10 km s“', which when taken at face 
value implies a mass ratio of Q = Mi/M 2 = 6 . 6 ^^ g (Icr), using the above value of X 2 . We derive constraints 
on the binary parameters from simultaneous modelling of the ellipsoidal light and radial velocity curves. We 
find 1(7 ranges for the photometric period (1.5430d < Rph < 1.5440d), X-velocity (350.2 < K 2 < 368.6 km s“'), 
inclination (67.0° < i < 77.4°), mass ratio (Q > 12.0), and orbital separation (1 1 . 55^0 < a < 12.5 QRq). Given 
these geometrical constraints we find the most likely value of the mass of the compact object is 9.41 Mq with a Icr 
range of 8.36M0 <M\< \0.16Mq. If we apply our tentative value of yrotsin/ = 90± 10 km s“* as an additional 
constraint in the ellipsoidal modelling, we find Icr ranges of 1.5432d < Rph < 1.5441 d for the photometric period, 
352.2 <K 2 < 370.1 km s“' for the X-velocity, 70.8° < i < 75.4° for the inclination, 6.1 < Q < ll.O for the mass 
ratio, and 12 . 35/^0 <a< 13.22 Rq for the orbital separation. These geometrical constraints imply the most likely 
value of the mass of the compact object of 1O.56M0 with a Icr range of 9 . 68 M 0 <Mi< 11.58M0. In either 
case the mass of the compact object is well above the maximum mass of a stable neutron star, and we therefore 
conclude XTE J1550-564 contains a black hole. 

Subject headings: binaries: spectroscopic — black hole physics — stars: individual (XTE J1550-564) — X-rays: 
stars 


1. INTRODUCTION 

X-ray novae provide the strongest evidence for the exis¬ 
tence of stellar mass black holes. These objects are interact¬ 
ing binaries containing a compact primary (a neutron star or 
a black hole) and what is usually a late-type secondary. X- 
ray novae spend most of their time in a low X-ray luminosity 
“quiescent” state, with an X-ray luminosity (Lft) roughly sim¬ 
ilar to the optical luminosity (Lopt), which includes a substan¬ 
tial contribution from the secondary star (Lstar)- The quiescent 


state is occasionally interrupted by “outbursts” where typically 
Lx Lopt ^ Tstar- In quiescence, the observed radial veloc¬ 
ity and light curves of the secondary star lead to dynamical 
mass estimates for the compact primary. If its mass exceeds 
the maximum stable mass of a neutron star (w 3Mq), the com¬ 
pact object is presumed to be a black hole (Chitre & Hartle 
1976; Kalogera & Baym 1996). In thirteen cases, the mass 
of the primary of an X-ray nova has been shown to exceed 
3M0, confirming the presence of black holes in these systems 


* Based on observations collected at the European Southern Observatory, Chile (program 67.D-0229), at the Magellan Walter Baade telescope at the Las Campanas 
Observatory, Chile, and the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica de Canarias 
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(Bailyn et al. 1995; Casares, Charles, & Naylor 1992; Casares, 
Charles, & Marsh 1995; Filippenko, Matheson, & Ho 1995; 
Filippenko et al. 1999; Filippenko & Chornock 2001; Greiner, 
Cuby, & McCaughrean 2001; McClintock & Remillard 1986; 
McClintock et al. 2001; Orosz et al. 1998, 2001; Remillard, 
McClintock, & Bailyn 1992; Remillard et al. 1996). These 
X-ray novae open up the possibility of studying the strong- 
held regime of general relativity. For example, quiescent X-ray 
spectra and ADAF (advection-dominated accretion how) mod¬ 
els provide evidence for black hole event horizons (Narayan, 
McClintock, & Yi 1996; Garcia et al. 2001). Also, the study of 
high-frequency QPOs (quasi-periodic oscillations) may lead to 
the hrst secure measurement of black hole spin (e.g. Remillard 
et al. 1999; Strohmayer 2001). 

XTF J1550-564 was discovered on 1998 September 7 by the 
All-Sky Monitor (ASM) aboard the Rossi X-ray Timing Ex¬ 
plorer (RXTE) (Smith et al. 1998). The optical counterpart (des¬ 
ignated V381 Normae) and the radio counterpart were discov¬ 
ered shortly thereafter (Orosz, Jain, & Bailyn 1998; Campbell- 
Wilson et al. 1998). This source was quickly identihed as a 
promising black hole candidate based on its rapid X-ray vari¬ 
ability, hard energy spectrum, and the absence of pulsations or 
X-ray bursts (Cui et al. 1999; Sobczak et al. 1999a). Relativis¬ 
tic plasma ejections at probable superluminal velocities were 
observed at radio wavelengths shortly after the strong X-ray 
hare in 1998 September (Hannikainen et al. 2001), indicating 
that XTF J1550-564 is another microquasar. The galactic mi¬ 
croquasars are excellent laboratories for the study of relativistic 
jets since they evolve orders of magnitude more quickly than do 
the jets in quasars (Mirabel & Rodriguez 1999). XTF J1550- 
564 is also of special interest owing to its complex X-ray vari¬ 
ability (e.g. Homan et al. 2001; Remillard et al. 2002). In this 
paper we report the results of our recent optical observations 
of XTF J1550-564. The observations and basic data reductions 
are summarized in §2. In §3 we establish the orbital parame¬ 
ters of the system, derive some properties of the secondary star, 
derive geometrical parameters for the binary, and discuss lim¬ 
its on the mass of the compact object. The implications of our 
results are discussed in §4 and summarized in §5. 

2. THE OBSERVATIONS AND THEIR REDUCTIONS 
2.1. Spectroscopy 

We obtained a total of 18 spectra of the source on 2001 May 
24-27 using Antu, which is the hrst 8.2 m telescope at the Eu¬ 
ropean Southern Observatory, Paranal. We used the FORSl 
imaging spectrograph with the 600V grism and a 077 wide slit; 
this combination gives a spectral resolution of 3.6 A FWHM 
and a wavelength coverage of 4571-6927 A. The seeing was 
better than 078 on May 24 and May 27, and relatively poor on 
May 26 (> 175). Clouds were present on the night of May 25, 
and only one useful spectrum was obtained in 175 seeing. The 
exposure times ranged from 20 to 40 minutes depending on the 
conditions. We also observed eight bright subgiant (luminos¬ 
ity class IV) and giant (luminosity class III) stars with spectral 
types from G8 to MO. An atmospheric dispersion corrector was 
used and the slit was kept at the default north-south direction. 
Following the standard procedure at Paranal, the hat-held and 
wavelength calibration exposures were obtained during the day¬ 
time hours with the telescope pointed at the zenith. 


We reduced all of the spectra using IRAF^. The standard 
tasks were used to apply the bias and Hat-held corrections, 
and to extract and wavelength-calibrate the spectra. The night- 
sky emission lines were used to make small adjustments to the 
wavelength scales of the spectra. The spectra were placed on 
an approximate Hux scale using observations (taken May 27) of 
the white dwarf hux standards EG 274 and Feige 110. The hnal 
reduced spectra of XTE J1550-564 have typical signal-to-noise 
ratios in the continuum near 6 per 1.16 A pixel at 6000 A and 
about 2 per 1.16 A pixel at 5000 A. 

We extended our grid of comparison stars by making use of 
archival observations. Eleven spectra of seven different stars 
were taken from the FORSl archive^. These spectra were ob¬ 
tained using the 600R grism and a 077 wide slit; this combi¬ 
nation gives a spectral resolution of 3.2 A FWHM and a wave¬ 
length coverage of 5142-7283 A. All of the standard reductions 
were performed with IRAF. We also used 49 spectra of 34 dif¬ 
ferent stars obtained 1995 April 30 and 1995 May 2-4 using 
the CTIO 4m telescope and the RC spectrograph and the Fo- 
ral 3072 x 1024 CCD (see Bailyn et al. 1995; Orosz & Bai¬ 
lyn 1997). The spectral resolution is 3.3 A FWHM and the 
wavelength coverage is 3850-7149 A. Finally, we used mod¬ 
erate resolution spectra of 6 stars with spectral types K5V to 
MOV obtained on 2001 February 16 and 17 with the red arm 
of the ISIS instrument on the 4.2m William Herschel Telescope 
(WHT) on Fa Palma. The combination of the 600R grating with 
the 1024 X 1024 TEK4 CCD and a 1 arcsecond slit width yields 
a resolution of w 70 km s“' FWHM and a wavelength coverage 
of 5856-6646 A. These spectra were reduced with IRAF. 

2.2. Photometry 

We observed XTE J1550-564 2001 June 26-28 using the 
Magellan Instant Camera (MagIC) on the Walter Baade 6.5m 
telescope located at Fas Campanas Observatory. The MagIC 
contains a SITe 2048 x 2048 CCD with a scale of 07069 per 
pixel. The source was observed for about 5 hours on each of 
the three nights in generally good conditions (0."3 to 0."7 see¬ 
ing and thin cirrus) with the Sloan r', and z' hlters. Some 
images in the Sloan g' filter were also obtained on the night of 
June 26. 

XTE J1550-564 has been monitored extensively by the Im 
YAFO telescope (Bailyn et al. 2000) at Cerro Tololo Interamer- 
ican Observatory since the discovery of the optical counterpart 
in 1998 September (Orosz, Jain, & Bailyn 1998). The V-band 
data from 2001 discussed here were collected and reduced in a 
similar manner to previous data from this source described by 
Jain et al. (2001b,c). 

As part of the normal acquisition procedure for spectroscopy 
we obtained a total of nine direct V images with exposure times 
of 30 to 90 seconds with VFT/FORS1 in its imaging mode. In 
addition, two additional direct images in B with exposure times 
of 120 and 200 seconds were taken. The image scale of FORSl 
is 072 per pixel, and the held of view is 6.8 x 6.8 square ar- 
cminutes. 

Finally, photometry of XTE J1550-564 was obtained 2001 
June 1 using the SuSI2 instrument on the Nasmyth focus of the 
3.5m New Technology Telescope (NTT) located at the Euro¬ 
pean Southern Observatory, Fa Silla. SuSI2 is a mosaic of two 


^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 

3 ESO programs 63.N-048I, 65.H-0360, and 265.D-50I6 
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2048 X 4096, thinned, anti-reflection coated EEV CCDs. We 
used 3x3 on-chip binning, yielding a scale of 0."24 per pixel 
and a held of view of 5.5 x 5.5 square arcminutes. The observ¬ 
ing conditions were quite poor: the seeing was between l."5 
and 2!'0 and there were passing clouds. We obtained a total 
of 51 usable /?-band images and 4 usable V -band images with 
exposure times of 3 to 5 minutes. 

The image processing routines in IRAE were used to cor¬ 
rect for the electronic bias and to apply the flat-field correc¬ 
tions to the VET, NTT, and Magellan images. The programs 
DAOPHOT IlE, ALLSTAR and DAOMASTER (Stetson 
1987; Stetson, Davis, & Crabtree 1991; Stetson 1992a,b) were 
used to compute the instrumental magnitudes of XTE J1550- 
564 and all of the field stars within about a 1.2 arcminute ra¬ 
dius. The comparison stars “A”, “B”, and “C” shown in Eigure 
2 of Jain et al. (2001b) were used to place the B- and V-band 
instrumental magnitudes on the standard scales. The R-band 
photometry and the Magellan photometry were left on the in¬ 
strumental magnitude scales. Since some of the V-band images 
from the VET and NTT were slightly under-exposed, we av¬ 
eraged groups of two to four consecutive exposures together, 
yielding two measurements from 24 May, and one measure¬ 
ment each from the nights of 25-27 May and 1 June. 

3. ANALYSIS 

3.1. Long-term Photometry 

Eigure la shows the complete YALO V-band light curve to¬ 
gether with our VET photometry. The 2-12 keV X-ray light 
curve from the All Sky Monitor (ASM) on RXTE is shown 
in Eigure lb. XTE J1550-564 had an optical reflare in 2001 
January (Jain, Bailyn, & Tomsick 2001a, Eig. la). Strangely 
enough, the corresponding outburst in X-rays was quite weak, 
with a peak level of less than 10 ASM counts per second. The 
source returned to its quiescent level by about April 20, and our 
V -band magnitudes from May 24-27 are fully consistent with 
the quiescent level: 21.83 < V < 22.24, with typical errors of 
0.04 mag. 

The source was not detected in B, and we place a conserva¬ 
tive lower limit of B > 24.0 ±0.1. Jain et al. (1999) had esti¬ 
mated B = 22.0 ± 0.5 for the quiescent level based on a detec¬ 
tion (near the plate limit) of the source on the SERC J survey 
print. The effective bandpass of the SERC J print is actually 
much redder than the standard Johnson B filter^, so the qui¬ 
escent magnitude given in Jain et al. (1999) is not a proper B 
magnitude but is closer to a V magnitude. 

3.2. Orbital Period and Spectroscopic Elements 

We measured the radial velocities of the secondary star us¬ 
ing th&fxcor task within IRAE, which is an implementation of 
the technique of Tonry & Davis (1979). The cross correlations 
were computed over the wavelength range 5000-6850 A, ex¬ 
cluding the Ha emission line, interstellar lines (Na D and the 
diffuse band near 5876 A), and a telluric feature near 6280 A. 
The spectra were continuum-subtracted prior to computing the 
cross correlation, and a Eourier filter was used to remove high 
frequency noise. With the exception of one spectrum from the 
end of the night of May 27 which had poor signal-to-noise, the 
cross correlation peaks were generally quite significant. The 
value of the Tonry & Davis “r” parameter, which is a measure 
of the signal-to-noise in a cross correlation, was generally in 

http://www.roe.ac.uk/ukstu/platelib.html#eqj 


the range of 3.0-4. 8 . The spectrum of the K3III star HD 181110 
usually gave the best cross correlation peaks as judged by the 
value of the r parameter, and the spectrum of the K4III star HD 
181480 generally gave the second-best cross correlation peaks. 
In the analysis below we adopt the velocities measured using 
the K3III template. Since the source was at or very near quies¬ 
cence, we believe our radial velocities are not biased by X-ray 
heating. 

To search for the spectroscopic period we computed a three- 
parameter sinusoidal fit to the 17 velocities for a dense range of 
trial periods between 0 and 4 days. The reduced values for 
these fits are shown in Eig. 2a. The free parameters at each trial 
period are the velocity semiamplitude K 2 , the epoch of maxi¬ 
mum velocity ro(spect), and the systemic velocity 7. The best 
fit is at a period of Pjp = 1.552 ±0.010 days (Icr error), where 
Xy = 1.38. The two alias periods near 0.7 and 2.5 days are ruled 
out by their large values of and by inspection of the folded 
velocity curves. We adopt the following spectroscopic elements 
with Icr errors: Psp= 1.552 ±0.010 days, X 2 = 349 ± 12 km s“', 
and ro(spect) = HJD 2,452,054.296 ± 0.014. In order to find 
the true systemic velocity we must know the radial velocity 
of the template star. We cross correlated the spectrum of HD 
181110 with 47 template spectra of stars with a radial veloc¬ 
ity measurement listed in the SIMBAD database. The median 
heliocentric velocity for HD 181110 was -71 km s“' and the 
standard deviation of the velocities was 18 km s“'. Based on 
this, we adopt 7 = -68 ±19 km s“' for XTE J1550-564. The 
velocities and the best fitting sinusoid are shown in Eig. 2b, and 
the spectroscopic elements are listed in Table 1. (The velocity 
plotted as an open circle in Eig. 2b was excluded from the fit 
since it deviates by more than 4 ct from the fit.) We note that our 
spectroscopic period is consistent with the photometric period 
of Pph = 1.541 ±0.009 days found by Jain et al. (2001b). The 
optical mass function is then 


f{M) = 




M\ sin^ i 


: 6.86 ± 0.71/W. 


©• 


( 1 ) 


27rG (Ml-1-/172)2 
(la error). The optical mass function places an immediate 
lower limit on the mass of the compact object, and since this 
minimum mass of the compact object is well above the max¬ 
imum mass of a stable neutron star (« 3Mq, Chitre & Hartle 
1976; Kalogera & Baym 1996), we conclude XTE J1550-564 
contains a black hole. See §3.5 for further discussion of the 
mass of the black hole. 


3.3. Parameters for the Secondary Star 

Using the spectroscopic ephemeris in Table 1, the eleven 
spectra with the highest signal-to-noise were Doppler corrected 
to zero velocity and averaged to create a “restframe” spectrum. 
The restframe spectrum was dereddened using the extinction 
law of Cardelli, Clayton, & Mathis (1989), where we adopted a 
visual extinction of Ay = 4.75 magnitudes [Nh = ( 8 . 5 ^ 24 ) x lO^' 
cm“2, Tomsick et al. (2001) and Ay =/Vh/1.79 x lO^', Predehl 
& Schmitt (1995)]. The true extinction is unlikely to be greater 
than 4.75 mag, since the extinction derived from X-ray mea¬ 
surements of Nh nearly always exceeds the optically measured 
value (Vtrilek et al. 1991). 

3.3.1. Spectral Type and Effective Temperature 

We used the technique outlined in Marsh, Robinson, & Wood 
(1994) to decompose the restframe spectrum into its disk and 
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stellar components. For this purpose we have three sets of tem¬ 
plate spectra: eight spectra of eight different stars taken with 
the FORSl 600V grism (resolution ~ 3.6 A), eleven spectra of 
seven different stars taken with the FORS1 600R grism (resolu¬ 
tion ~ 3.6 A), and 49 spectra of 34 different stars taken at CTIO 
(resolution ~ 3.3 A). The dereddened restframe spectrum and 
the template spectra were all normalized to unity at 6000 A and 
resampled to a pixel size of 1.155 A. The template spectra were 
then scaled by various values of a weight factor w between 0.0 
and 1.0 in steps of 0.02, and subtracted from the restframe spec¬ 
trum. The scatter in each difference spectrum was measured by 
fitting a low-order polynomial and computing the rms differ¬ 
ence over the same wavelength region used for the cross cor¬ 
relation analysis above. For each template spectrum we looked 
for the value of w that gave the “smoothest” difference spectrum 
(i.e. the w with the lowest associated rms). Figure 3 shows the 
rms values as a function of the template spectral type for the 
three sets of templates. Since the spectral resolution in each of 
the three sets is slightly different, one should only compare the 
rms values within a single set. For the CTIO grid the minimum 
rms is near spectral type G6. However, the difference spectra 
for the G6 templates (and the templates of earlier type) have 
negative flux in the blue parts of the spectrum, which is un¬ 
physical. This difficulty could be lessened if one were to adopt 
a larger value of Ay; however, values larger than Ay = 4.75 mag 
are unlikely (see above). The templates with spectral type G8 
and later, on the other hand, all produce subtracted spectra with 
positive differences. The minimum rms values for the FORSl 
templates is at spectral type G8 for the grism 600V and spectral 
type K1 for the grism 600R. In general, there is a steep rise in 
the rms values for template spectral types later than K5 and ear¬ 
lier than G5. The weight factor w is roughly correlated with the 
template spectral type where templates with later spectral types 
have smaller values of w. 

Fig. 4 shows the restframe spectrum decomposed into the 
disk and stellar components for two different templates ob¬ 
served with the FORSl grism 600V (i.e. the same instrumental 
configuration used for the XTF J1550-564 spectra). Fig. 4a 
shows the decomposition using the K3III template HD 181110 
(w = 0.48) and Fig. 4b shows the decomposition using the G8IV 
template HD 157931 (w = 0.74). In both cases the disk spec¬ 
trum is basically flat in F\, and the only obvious features re¬ 
maining are the Ha emission line, the interstellar lines, and the 
telluric line (note that these features have been smeared in the 
process of making the restframe spectrum). The Mg b lines 
near 5169 A are matched better with the G8IV template. Apart 
from Mg b, the two difference spectra look very similar to the 
eye. 

Based on the decomposition analysis above, we adopt a spec¬ 
tral type of G8-K4 for the secondary star. The spectral type is 
almost certainly not K7 or later since the TiO bands near 6230 
and 6700 A become apparent in dwarfs at K7 and later and in 
giants at K5 and later, and these bands are not evident in our 
spectrum. Spectral types earlier than about G5 are also firmly 
ruled out based on the high rms values (Fig. 4b) and the partly 
negative difference spectra. The secondary star cannot have a 
normal main sequence gravity owing to the relatively long or¬ 
bital period, so its surface gravity must be somewhat less than 
the nominal values for luminosity class V. The surface gravity 
of the secondary is easily found using Kepler’s Third Faw, as¬ 
suming that the star fills its Roche lobe. In that case, there is 
a relationship between the mean density of the secondary star 


and the orbital period: 

Ml Stt 1 

^ IttRI P^GRUQ)1 + Q 

where Rrl{Q) is the sphere-equivalent radius of the Roche lobe 
for unit separation, and G is the universal constant of grav¬ 
ity. One can readily show through numerical integration that 
for 2 < Q < 20, the quantity R^i(Q)(l + QT^ only varies by 
4.6% (between 9.7 and 10.2). Hence, for a fixed orbital pe¬ 
riod P, the mean density of the secondary is nearly indepen¬ 
dent of the mass ratio Q. This in turn means that the result¬ 
ing surface gravity is a very weak function of the assumed 
secondary star mass Mi (logg oc 5 logM 2 ). We find surface 
gravities of logg = 3.51, 3.61, and 3.72 for secondary star 
masses of 0.5, 1.0, and 2.0Mq, respectively. These values of 
logg imply a luminosity class between III and V [e.g. a K4V 
star nominally has logg = 4.60 and a K4III star nominally has 
logg = 1.7, Gray (1992)]. The temperature corresponding to 
a given spectral type depends on the luminosity class, where 
main sequence stars of a given spectral type (i.e. luminosity 
class V) are hotter than the corresponding giants (luminosity 
class III). Straizys & Kuriliene (1981) give Teff ~ 4100 K for a 
K4III star and Teff ^ 5100 K for a G8IV star. Thus we adopt 
4100 < Teff < 5100 K for the analysis below. 

3.3.2. Rotational Broadening 

We used the Marsh et al. (1994) technique to estimate the 
projected rotational velocity Wot sin/ of the XTF J1550-564 
restframe spectrum. The dereddened restframe spectrum, the 
spectrum of the G8IV star HD 157931, and the spectrum of the 
K3III star HD 181110 were each normalized to unity at 6000 A 
and resampled to a pixel size of 0.46 A. An IRAF procedure im¬ 
plementing the analytic broadening kernel given in Gray (1992, 
p. 374) was used to broaden the template spectrum using values 
of Wot sin / between 50 and 160 km s“' in steps of 5 km s“* (we 
adopted a linear limb darkening coefficient of 0.6). The value 
of the weight factor w was optimized to find the minimum rms 
value at each input value of Wot sin/. The results are shown in 
Fig. 5a. The minimum rms value occurs for Wot sin / = 90 km 
s“' for the G8IV template and at Wot sin / = 95 km s“' for the 
K3III template. 

We next attempted to measure Wot sin / for the secondary star 
using a different technique that is based on the relationship be¬ 
tween the width of the cross correlation (cc) peak and the ro¬ 
tational velocities of the object and template spectra (Horne, 
Wade, & Szkody 1986; Fchevarria et al. 1989; Casares et al. 
1993). A Gaussian with a fitting width of 750 km s“' and 
a baseline of zero was used to determine the width of the cc 
peaks. The measured widths were relatively insensitive to the 
exact width of the fitting region. The full width at half maxi¬ 
mum (FWHM) of the restframe spectrum cc peak is in the range 
of 285 to 305 km s“', depending on which template is used. In 
order to find out what rotational velocity this corresponds to, 
we constructed calibration curves using some bright compari¬ 
son stars observed with the FORS 1 600V grism. Three differ¬ 
ent comparison stars were each broadened by various amounts, 
and the FWHM of the cc peaks were measured using the K3III 
template HD 181110. The results are shown in Fig. 5b. These 
calibration curves show that a cc peak FWHM of 305 km s“' 
corresponds to a rotational velocity of Wot sin / w 90 km s“', and 
that a FWHM of 285 km s“' corresponds to an input Wot sin / of 
about 50 km s“'. 
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The Marsh et al. (1994) technique indicates a rotational ve¬ 
locity of Vrotsin/ = 90 km s“' (given the better performance 
of the G8 template in separating the star and disk spectra we 
adopted the value of Vrot sin i measured with this template, Fig. 
5a), and the cross correlation technique indicates Vrotsin/ < 90 
km s“' (Fig. 5b). Based on the results of the two techniques, 
we adopt a value of Vrot sin / = 90 ± 10 km s“' for the rotational 
velocity of the secondary star. 

The value of Vrotsin/ should be treated with caution since 
our spectral resolution is relatively low (about 160 km s“') and 
the signal-to-noise ratio is modest (about 20 per pixel in the 
restframe spectrum). In order to have some idea of the size 
of the systematic errors involved with measuring the rotational 
velocity in a relatively noisy low resolution spectrum, we per¬ 
formed some numerical experiments using the moderate reso¬ 
lution spectra obtained with the WHT. These experiments are 
by no means comprehensive since we have only a limited range 
of spectral types (K5 to MO, all dwarfs) and wavelength cov¬ 
erage (roughly between the Na D lines and Ha). The six tem¬ 
plate spectra were first normalized at a common wavelength 
(6450 A). A simulated object spectrum was made as follows: 
a template was scaled by 0.7, rotationally broadened by 90 km 
s“\ degraded to a resolution of 3.6 A by convolution with a 
Gaussian, and resampled to 1.16 A (to match the pixel size of 
the VLT/FORSl spectra). Next, a constant value of 0.3 was 
added (to simulate the accretion disk), and Gaussian noise was 
added so that the signal-to-noise ratio was « 20. Finally, the 
simulated object spectrum and the remaining template spectra 
were resampled to a pixel size of 0.46 A. The rotational ve¬ 
locity of the simulated object spectrum was extracted by the 
two different techniques, and the results are shown in Figs. 5c 
and d. The Marsh et al. technique gives a rotational velocity 
of just over 100 km s“' using the two best matched templates, 
and the cross correlation technique gives a rotational velocity 
between about 80 and 110 km s“' using all templates and be¬ 
tween about 100 and 110 km s“' using the two best matched 
templates. Given that the input rotational velocity for the sim¬ 
ulated spectrum was 90 km s“', we conclude that the system¬ 
atic errors involved in this technique are not excessively large 
and are probably on the order of our adopted Icr error of 10 
km s“'. Needless to say it would be desirable to repeat the 
measurement of the rotational velocity of the secondary star in 
XTE J1550-564 using data with better resolution and a higher 
signal-to-noise ratio. 

If the secondary star is rotating synchronously with the orbit, 
its mean rotational velocity depends only on the binary mass 
ratio Q = M1/M2 and the projected orbital velocity K2 (e.g. 
Horne, Wade, & Szkody 1986). We used the ETC code (Orosz 
& Hauschildt 2000) to numerically compute sphere-equivalent 
Roche lobe radii as a function of the mass ratio Q, and a simple 
Monte Carlo code to propagate the (Icr) uncertainties on the 
input parameters. We find Q = 6.6j;[ g (Ict) for K 2 = 349 ± 12 
km s“' and Vrotsin/ = 90± 10 km s“'. This mass ratio is rela¬ 
tively far from unity, which is typical of X-ray novae with cool 
low mass companions (e.g. Bailyn et al. 1998 and references 
therein). 

3.4. The Ha Profile 

The XTE J1550-564 spectra show a strong, double-peaked 
Ha emission line. Eig. 6 shows the mean Ha line profile ob¬ 
tained by averaging the eleven spectra with the highest signal- 
to-noise. The equivalent width of the line is about 24.7 A, the 


full width at half maximum is about 1500 km s“', and the max¬ 
imum velocity in the line wings is about 1000 km s“'. We do 
not detect any feature at the wavelength of H/3, although the 
signal-to-noise is rather poor there. 

The double-peaked emission line is characteristic of emis¬ 
sion from an accretion disk (Smak 1981). In this case there 
is < 8% stellar contribution to the Ha profile since late G and 
early K stars have relatively weak Ha absorption lines (EW 
< 2 A, Eig. 4). We modelled the line profile using a simple 
model outlined in Smak (1981; see also Johnston, Kulkarni, & 
Oke 1989, and Orosz et al. 1994). The three main parameters 
are (1) the power-law exponent a, which determines the density 
function of the emitting atoms, /(r) oc r““ (where the coordi¬ 
nate r is normalized to 1 at the outer edge of the disk), (2) the 
ratio ri of the inner disk radius to the outer disk radius, and (3) 
the radial velocity Vd of the outer edge of the disk. Horne & 
Marsh (1986) have modified the model to include profiles from 
an optically thick disk; this modification requires that the in¬ 
clination be included as an additional parameter. An optically 
thin model with a = 2, r\ = 0.17, and Vd = 420 km s“' pro¬ 
vides a reasonably good match to the observed profile (we have 
made no attempt to optimize the model fits), although there is 
some excess emission in the blue line wing and the central de¬ 
pression in the observed profile is deeper than the model (see 
Eig. 6). The optically thick models have deeper central depres¬ 
sions, and the depth of an optically thick model with i = 50° 
(and a = 2, ri = 0.17, and Vd = 420 km s“', Eig. 6) matches the 
depth of the central depression reasonably well. One should 
note that the inclination used in the model profile may not be a 
good indicator of the actual binary inclination (see §3.5). Eor 
example, the depth of the central depression in the Ha profile of 
GS 1124-683 (XN Mus91) was observed to vary over an orbital 
cycle (Orosz et al. 1994). 

The line wings of the Ha profile of XTE J1550-564 are rel¬ 
atively steep and do not extend to very high velocities (slightly 
more than 1000 km s“'). In contrast, the line wings of the Ha 
profile of GS 1124-683 and H1705-250 (XN Oph77) extend 
out to almost 2000 km s“* (Orosz et al. 1994; Remillard et al. 
1996). Eor A0620-00, the width of the Ha profile has been 
observed to vary between these extremes during 1991-1993 
(Orosz et al. 1994). The weakness or absence of the H/3 line 
in XTE J1550-564 is similar to what is observed in H1705-250 
in 1994 (Remillard et al. 1996). There is some excess emission 
in the blue line wing of the XTE J1550-564 Ha profile. Excess 
emission in the blue line wing was observed in GS 1124-683 
in 1992, but not in 1993, and in A0620-00 in 1991, but not in 
1993. Einally, we note that the ratio of the disk velocity Vd to 
the X-velocity in XTE J1550-564 is = 1.2. This value is 
well within the range observed for the other black hole X-ray 
novae with (cool) low mass companions (see Table 2). Thus 
in terms of the properties of its quiescent accretion disk, XTE 
J1550-564 appears to be a typical black hole system. 

3.5. Ellipsoidal Light Curve Models and Derived 
Astrophysical Parameters 

The value of the orbital inclination i is needed to compute the 
component masses and other interesting parameters. Normally 
one models the ellipsoidal light curves to estimate the inclina¬ 
tion of a semi-detached system such as XTE J1550-564, and 
for this purpose we used the EEC code (Orosz & Hauschildt 
2000). The model has many free parameters. The binary pa¬ 
rameters for a circular orbit are the inclination i, the mass ratio 
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Q, the orbital separation a, the orbital period P, and the epoch 
of the inferior conjunction of the secondary star 7o(photo). The 
main parameters for the secondary star are its mean temperature 
Teff, its gravity darkening exponent /?, its Roche lobe filling fac¬ 
tor /, and its rotational velocity with respect to synchronous O. 
Finally there are parameters related to the accretion disk which 
are the outer radius of the disk rout, the inner radius of the disk 
Tin, the opening angle of the outer disk rim the temperature 

at the inner edge Tdisk, and the power-law exponent on the disk 
temperature profile where T{r) = 7disk(r/rin)^. Some of these 
parameters can be fixed at reasonable values. We assume the 
secondary star exactly fills its Roche lobe and is in synchronous 
rotation. Hence f = fl= 1. Following Claret (2000) the value of 
the gravity darkening exponent /3 is a function of the mean tem¬ 
perature of the secondary star. In this case the gravity darken¬ 
ing exponent is in the range of 0.06 < /3 < 0.11 for our adopted 
temperature range of 4100 < Teff < 5100 K. We used the spe¬ 
cific intensities from the NextGen models (Hauschildt, Al¬ 
lard, & Baron 1999a; Hauschildt et al. 1999b), so consequently 
no parameterized limb darkening law was needed. The inner 
radius of the disk was taken to be 3% of the compact object’s 
Roche lobe radius, which is roughly lO'* Schwartzshild radii 
(see Hameury et al. 1997). 

We are left with ten free parameters in the ellipsoidal model: 
Q, i, a, P, 7’o(photo), Teff, Tdisk, ?-out, and (3n^. For a given 
ellipsoidal model, we can compute various observable quanti¬ 
ties such as the shape of the light curves in a given bandpass, 
the amplitude and phasing of the radial velocity curve of the 
secondary star, the amount of disk light in various bandpasses, 
the radius, gravity, and rotational velocity of the secondary star, 
etc. The goal is to find the set of model parameters which gives 
the best fit to all the observed properties of XTF J1550-564. 

To accomplish the optimization goal we used a genetic op¬ 
timization code based on the PIKAIA routine given in Char- 
bonneau (1995), with the “black sheep” modification outlined 
by Bobinger (2000). This genetic code is very efficient at find¬ 
ing the global optimal solution in a large parameter space (see 
Charbonneau 1995 and Metcalfe 2001 for more detailed discus¬ 
sions of genetic algorithms). The algorithm used in the code is 
quite simple. We define a “population” of 100 random param¬ 
eter sets. For each parameter set we compute the ellipsoidal 
model and define a “fitness” based on how well the particular 
model compares with the available observations. After the fit¬ 
ness of the initial population is determined, 100 new parameter 
sets are produced as a result of “breeding” between pairs of 
members, where the probability of breeding is based on the fit¬ 
ness. Random variations (i.e. “mutations”) are introduced into a 
small fraction of the breeding events. The process of breeding a 
new population and evaluating its members is repeated (i.e. the 
population “evolves”) over many generations (typically a few 
hundred or more) until convergence is achieved. 

We will use the statistic to evaluate the goodness-of-fit 
between the predicted and observed light and velocity curves: 

2 _ 2 , 2 
Xtotal “ Xlight Xvelocity 




(y(x,;ai...flio)-y,)^ 


. (y(x,;ai...aio)-y;)^ 


Here the notation y(i,;ai.. .aw) means the “model value com¬ 
puted at X,” and y, means the “observed quantity at the same 
Xi.” There is a double summation in the Xught because we 
have six filters. In the genetic code, the fitness of a population 
member is related to its relative value of the x^- The member 


with the lowest x^ has the highest fitness, the one with the sec¬ 
ond lowest x^ has the second highest fitness, and so on. We 
can also consider at least three additional measures of fitness: 
(i) For each model we can predict 0, the duration of an X-ray 
eclipse, if any. In the case of XTF J1550-564, the X-ray eclipse 
duration is 0 = 0 at a high level of confidence since there were 
no eclipses observed in the almost 300 pointed RXTE observa¬ 
tions of XTF J1550-564 during its various outbursts nor in the 
extensive RXTE ASM observations (Fig. lb). Thus we have: 


2 _f0, if 0 = 0; 
^e“ll0®, if0>O. 


(3) 


As a result of this definition of Xs, any model where the X-ray 
source is eclipsed is assigned a very low fitness, (ii) For each 
model we can predict ky, which is the fraction of the V-band 
light contributed by the accretion disk. Observationally, ky can 
be estimated using the decompositions discussed in §3.3. At 
5500 A, the disk fraction is w 0.4 for the K3III template (Fig. 
4a), and « 0.2 for the G8 IV template (Fig. 4b). For the pur¬ 
poses of the ellipsoidal modelling we will adopt ky = 0.3 ± 0.1. 
We then have 


2 _ (ky(ai.. .aio)-0.3)^ 
X<:y ~ 0 1^ 


For simplicity we will neglect the small correlation between the 
disk fraction and the effective temperature, (iii) For each model 
we can predict Vjot sin i, which is the value of the projected ro¬ 
tational velocity of the secondary star. In §3.3.2 we tentatively 
measured Vot sin / = 90 ± 10 km s“'. We then have 


2 _ (Trotsin i(ai ...aw)- 90)^ 

Xrot - 2 q2 


Our total becomes xLai = Xught + xL + Xe + xl + X?of 

In any optimization procedure, one must decide how to as¬ 
sign weights to the observations. In our case some initial fits 
were performed to get an approximate solution. Based on these 
initial results, the error bars on the photometric measurements 
were scaled slightly to give Xught/ ~ 1 at the minimum in each 
of the six filters separately. Similarly, the errors on the radial 
velocities were scaled to give Xrei/^ w 1 at the minimum. 

Since our measured value of the rotational velocity may have 
a large systemic error (the value we measured is smaller than 
the spectral resolution), we ran two sets of fits. Fitting run A 
did not include the contribution of Vrotsin/ in the total x^- Fit¬ 
ting run B did include the observed value of Wot sin i = 90 km 
s“^ in the total x^ (i-C- we are taking the value of 90 km s“' at 
face value). For each of the two fitting runs A and B, the ge¬ 
netic fitting code was run six separate times on the data with the 
scaled error bars, using a population size of 100. The volume 
of parameter space that was searched was the same in all six 
cases, but the order in which the parameters were encoded into 
the “genes” was different. In practice this means that the initial 
random populations were different, and the way two population 
members exchange parameters when breeding was different. In 
all cases the final best-fitting parameters were statistically iden¬ 
tical after about 200 generations, giving us great confidence that 
our solutions for each of runs A and B are truly the global ones. 

The folded light curves in the six bands and the model light 
curves are shown in Fig. 7. Although we fit the unbinned data, 
we also show in Fig. 7 smoothed versions of the light curves 
which were made by computing the magnitude within bins 0.05 
phase units wide. The smoothed light curves generally follow 
the fitted models quite well. 
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In order to estimate the uncertainties on the fitted parameters 
and the uncertainties on the derived astrophysical parameters 
we collapse the 10-dimensional function onto the appropri¬ 
ate parameter of interest. The 1, 2, and 3 (t confidence limits are 
then the values of the parameter where = x^i„ +1 , Xmin 
and Xmin+ 9, respectively. Using the genetic fitting code and 
the Black Sheep modification it is quite easy to compute these 
confidence limits since the region of parameter space near the 
X^ minimum is sampled quite well. For each model that is 
computed, the x^ of the fit, the values of the free parameters, 
and the computed astrophysical parameters (e.g. the black hole 
mass, the secondary star radius, etc.) are saved. After a suffi¬ 
ciently large number of models have been computed (more than 
515,000 in fitting run A and 465,000 in fitting run B), we can 
collapse the x^ function onto the parameter of interest by sim¬ 
ply plotting the x^ vs. the parameter of interest and looking at 
the “lower envelope” (i.e. the minimum x^ within small bins 
over the whole range). 

Fig. 8 shows the collapsed x^ function for the ten fitted pa¬ 
rameters in the fitting run A, and Fig. 9 shows the similar plot 
for the fitting run B. In some cases, we fitted fifth or sixth order 
polynomials to the curves to more objectively to determine the 
value of the parameter where x^ = Xmin '^he 1, 2, and 3 a 
confidence limits. In general, the important geometrical param¬ 
eters in the ellipsoidal model (i.e. the inclination i, separation a, 
etc.) have well-defined minima and 3a limits. In addition, the 
derived parameter ranges are not too different between fitting 
run A and fitting run B. The one exception is the mass ratio Q. 
In the fitting run A, where the observed rotational velocity was 
not included in the total x^, the value of Q can only be assigned 
a lower limit (Q > 6 at the 3a level). However, the mass ratio Q 
does have well defined upper and lower limits in the fitting run 
B, where the observed rotational velocity was included in the 
total x^- On the other hand, the accretion disk parameters and 
the effective temperature of the secondary star are not so well 
constrained. 

Fig. 10 shows the collapsed x^ function for six derived pa¬ 
rameters in fitting run A; Mi (the compact object’s mass), and 
M 2 (the secondary star’s mass), Uotsin/, (the predicted rota¬ 
tional velocity of the secondary star), R 2 (secondary star’s ra¬ 
dius), logg (the secondary star surface gravity), and La (the sec¬ 
ondary star’s bolometric luminosity). Fig. 11 shows the similar 
plot for fitting run B. Key results are summarized are in Table 
3. The nominal mass of the black hole is 9.41 Mq for fitting run 
A and 1O.56M0 for fitting run B, respectively. We will adopt 
the results from fitting run B in the following discussions. 

The only observational constraint related to the disk at 
present is that the disk should contribute roughly 30% of the 
light in the V-band, and there is a wide range of accretion 
disk parameters which produce a disk with roughly the correct 
brightness. Since the inclination seems to be fairly high, the 
disk may partially eclipse the secondary star, and vice versa. 
Unfortunately, we do not have complete phase coverage, so 
we cannot at present detect the subtle signatures of grazing 
eclipses. Therefore, observations with much better phase cov¬ 
erage near the conjunction phases should be obtained. Al¬ 
ternatively, observations in the infrared (i.e. J-, H-, and K- 
bands) would also be especially useful since the usual assump¬ 
tion is that the relative contribution of the accretion disk is much 
smaller at these wavelengths. This assumption can and should 
be tested with good quality infrared spectra. If the disk turns out 

^ ftp://calvin.physast.uga.edu/pub/NextGen/Colors/ 


to be faint in the infrared, then the model fits would of course 
be much less sensitive to the assumed disk parameters. 

The inclusion of the rotational velocity as an additional fit¬ 
ting constraint resulted in only a modest (« 10%) change in 
the derived mass of the black hole. This is because the black 
hole mass scales as /(M)(l + l/Qf, and the mass ratio Q is 
relatively large. On the other hand, the derived properties of 
the secondary star are quite different between the two fitting 
runs. For example, M 2 < 0.79Mq (la) for fitting run A and 
M 2 = I. 3 IM 0 for fitting run B. Clearly, a better measurement 
of Wot sin i would be quite useful since it would place a much 
stronger constraint on the mass ratio, which in turn would place 
stronger constraints on the properties of the secondary star. In 
order to significantly improve the precision of the black hole 
mass, one would need to also improve the precision of the mass 
function f(M). 

Finally, it is worthwhile to point out that the results sum¬ 
marized in Table 3 are our best estimates based on the data 
we currently have. Although the parameter uncertainties have 
been rigorously computed, they represent statistical errors only. 
There may be systematic errors in, for example, the measured 
rotational velocity Uot sin / owing to the low spectral resolution 
or in the derived inclination owing to the less than perfect phase 
coverage and errors in measuring the contamination from the 
accretion disk. It was straightforward to perform two fitting 
sequences to see what effect including the rotational velocity 
as a fitting constraint had on the outcome. It is less straight¬ 
forward to quantify the effects of incomplete phase coverage 
or imperfect disk/stellar decomposition. However, we can say 
with some confidence that since our derived inclination is al¬ 
ready quite high and near the upper limit imposed by the lack 
of X-ray eclipses, it is unlikely that we have underestimated the 
true inclination significantly. Ultimately one would need to ob¬ 
tain additional data to see if similar light curves are observed. 

3.6. Distance 

The distance to the source is a function of four parameters: 
The mass M 2 and temperature Teff of the secondary star fix its 
luminosity (see Eq. 2); interstellar extinction (Ay) makes the 
source appear fainter; and extra light from the accretion disk 
(parameterized by the V-band “disk fraction” ky) makes the 
source appear brighter. We have reasonably good values for 
all of these parameters. As noted in §3.3, Ay « 4.75 can be 
derived by modelling the X-ray spectrum and assuming a typi¬ 
cal gas to dust ratio for the interstellar medium; ky « 0.30 was 
determined from our spectral decomposition (Figs. 4a and 4b); 
and a temperature range of 4100 < Teff < 5100 K was deduced 
for the secondary star. Using these results and a simple Monte 
Carlo procedure, we computed the distance and its uncertainty 
for various values of the secondary star mass /W 2 . We used the 
synthetic photometry computed from the NextGen models^ to 
determine the expected absolute V magnitude of the star from 
its temperature, radius, and surface gravity. For this procedure 
we adopted 4100 < Teff < 5100 K, ky = 0.30 ±0.05, a mean ap¬ 
parent V magnitude of 22.0 ± 0.2, and Nh = (8.5 ± 2.3) x 10^^ 
cm“^ (90% confidence) and Ay =Nh/1-79 x 10^^ (exact). The 
results are displayed in Table 4 for a wide range of secondary 
star masses. For masses in the range of 0.15 to 3.QMq the com¬ 
puted distance is in the range 3.0 to 7.6 kpc. However, at a fixed 
mass, the formal error on the distance is rather large due to the 
relatively large adopted temperature range and to the relatively 
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large error in Nh- The distance for the case of M 2 = 1.31 Mq is 
5.9 kpc, again with a rather large formal error (Icr range of 2.8 
to 7.6 kpc). 

We can make an independent estimate of the distance if we 
assume the observed systemic velocity of 7 = -68 ±12 km s“^ 
is entirely due to differential galactic rotation. We use the rota¬ 
tion curve of Fich, Blitz, & Stark (1989) and the standard lAU 
rotation constants of /?o = 8.5 kpc and 0o = 220 km s“'. The ro¬ 
tation curve for a source with galactic coordinates of f-325.88° 
and b = -1.83° is roughly parabolic, with a minimum radial ve¬ 
locity of « -95 km s“' at a distance of w 7 kpc. The expected 
distance for a source with 7 = -68 km s“^ is either w 4.4 kpc 
or w 9.7 kpc. Considering the errors on 7, the kinematic dis¬ 
tance can be anywhere between « 3.2 kpc and « 10.8 kpc at the 
1(7 level, which is entirely consistent with the ranges derived 
above. Given the general agreement between the two distance 
estimates we conclude that XTE J1550-564 has essentially no 
(radial) velocity with respect to its local standard of rest. 

We can use the distance estimates for XTE J1550-564 to 
calculate an order of magnitude estimate for the maximum X- 
ray luminosity seen during outbursts. The spectral analysis of 
Sobczak et al. (2000) for the intense flare of 1998 September 
19 implies an unabsorbed flux of 3.7 x 10“^ erg cm“^ s“^ over 
the range of 1-20 keV, with the spectrum dominated by the 
power-law component. For our favored distance of 5.3 kpc, the 
isotropic luminosity is then 1.77 x 10^® erg s“^ (with a large 
error owing to the uncertainty in the distance). This is just 
slightly above the Eddington luminosity for a 1O.56M0 black 
hole, which is a remarkable coincidence given the significant 
uncertainty in the distance. 

4. DISCUSSION 

Bailyn et al. (1998) analyzed the distribution of black hole 
masses in seven systems and found a high probability that six of 
the seven systems (GRO J0422H-32, A0620-00, GS 1124-683, 
GRO J1655-40, H1705-250, and GS 2000 h- 25) have masses 
which are consistent with w 1 Mq. The seventh system, V404 
Cyg, has a mass which is significantly larger (about I 2 M 0 ; 
Shahbazet al. 1994). The mass of the black hole inXTEJ1550- 
564 is well above 1 Mq. With My > 8 .IM 0 at 3 (t confidence 
(Table 3), its mass may in fact be similar to the mass of the 
black hole in V404 Cyg. Since the work of Bailyn et al. (1998), 
the mass functions for seven additional systems have been mea¬ 
sured, and improved parameters for some of the original seven 
systems have been measured. Thus the issue of the observed 
black hole mass distribution should be revisited to see if the 
clustering of black hole masses near 7 Mq is still significant. 
Given the new data, we can begin to make meaningful compar¬ 
isons with formation theory. For example, the detailed forma¬ 
tion models of Fryer & Kalogera (2001) predict a mass distribu¬ 
tion which is continuous and which extends over a broad range 
(in particular, they predict no peak at « 1 Mq). 

The determination of the black hole mass for XTE J1550- 
564 is especially important for the interpretation of the high- 
frequency X-ray quasi-periodic oscillations (QPOs) observed 
for this system (Remillard et al. 2002; Homan et al. 2001). 
Models for several types of oscillations predicted in general rel¬ 
ativity are under investigation as possible causes of these QPOs 
(e.g. Remillard 2001, and references therein); all of these de¬ 
pend on both the mass and spin of the black hole, and possibly 
also on conditions in the inner accretion disk. Despite con¬ 
siderable uncertainties in the models, we can offer a few com¬ 


ments on the implications of our mass determination for XTE 
J1550-564. At the nominal mass of 1 O. 7 M 0 , oscillations at 
the frequency of the last stable orbit for a Schwarzschild black 
hole (spin parameter, < 2 * = 0 ) would be seen at a frequency, 
= 2199 /(Mbh/M 0) = 208.2 Hz (Shapiro & Teukolsky 1983). 
Since XTE J1550-564 has exhibited QPOs with frequencies up 
to 284 Hz (Homan et al. 2001), it appears that plausible mech¬ 
anisms require a* > 0. (The oscillation frequency would be 
284 Hz (for a* = 0) if Mbh = 1 .IAMq, which is below the 3 (j 
lower limit on the black hole mass.) This argument was used 
by Strohmayer (2001) to suggest that the 450 Hz QPO in GRO 
J1655-40 implies appreciable spin for the black hole in that sys¬ 
tem. We further note that the ratio of the black hole masses for 
XTE J1550-564 (Mi = lO.bMg, Table 3) and GRO J1655-40 
(Ml = 6.3Mq, Greene, Bailyn, & Orosz 2001) is ~ 1.7, which 
is nearly inversely proportional to their maximum QPO fre¬ 
quencies: (284Hz/450Hz)“^ ~ 1. 6 . We may then speculate that 
these results are consistent with a common QPO origin, with 
QPO frequencies that vary as M“\ which might be expected if 
the black holes have similar values of the spin parameter. 

Both XTE J1550-564 and GRO J1655-40 have had several 
closely-spaced outburst events relatively shortly after their ini¬ 
tial discoveries. Interestingly enough, optical precursors to X- 
ray outbursts (i.e. X-ray delays) have been observed for both 
sources. According to Hameury et al. (1997), these X-ray de¬ 
lays can best be understood using a two component accretion 
model consisting of a standard thin disk in the outer regions and 
an ADAF region in the inner region. At the onset of an outburst, 
the heating front reaches the ADAF region relatively quickly, 
but then must propagate more slowly on a viscous timescale, 
thereby greatly delaying the production of X-rays. Orosz et al. 
(1997) observed an optical precursor to the 1996 April outburst 
of GRO J1655-40. In that case, the delay between the V-band 
rise and the detection by the ASM on RXTE was 5.6 ±0.8 days. 
In the case of the 2000 outburst of XTE J1550-564, Jain et al. 
(2001c) reported a delay of 8.8 ±0.6 days between the rise in V 
and the detection by the RXTE ASM. The ratio of these delays 
is 8 .8/5 .6 ~ 1.6, which is nearly the same as the ratio of the 
masses and is the same as the inverse ratio of their maximum 
QPO frequencies (~ 1.7 and ~ 1 .6 respectively). We may again 
speculate that these results imply that the size of the ADAF cav¬ 
ity scales with the mass of the black hole. 

For the case of a companion star mass of \.3Mq (Table 3), 
the companion star in XTE J1550-564 appears to have mass 
that is larger than the typical mass for X-ray novae with a cool 
(K or M) companion. For comparison, the KO secondary star 
in V404 Cyg has a mass of about 0.1 Mq (Shahbaz et al. 1994). 
Since the secondary star in XTE J1550-564 is about three times 
larger than a normal K3/4 dwarf, its evolutionary state may be 
different than the dwarf secondaries in the short period X-ray 
novae such as A0620-00 or GS 1124-683. If the star is still 
in the core hydrogen burning phase, then its nuclear evolution 
timescale would be « 1 x 10 ^ years, leading to an average mass 
transfer rate of 10“^ M 0 yr“/ If the star has finished its core 
hydrogen burning, then presumably its radius is expanding on 
a thermal timescale of w 10 ^ years, leading to an average mass 
transferrate of \0~''Mq yr“^ 

We can make a rough estimate of the mass transfer rate in 
XTE J1550-564 as follows. The RXTE satellite performed nu¬ 
merous observations of XTE J1550-564 during its 1998 and 
1999 outbursts. Sobczak et al. (2000) give the unabsorbed 2- 
20 keV flux and the 20-100 keV flux for the « 200 PC A and 
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HEXTE observations from this period, respectively (their Eig. 
6 ). Integrating under the curves, we find a total 2-100 keV flu- 
ence of 1.06 erg cm“^ for the 1998-1999 outburst. Integrat¬ 
ing the ASM light curve (Eig. 1) we find that the outburst in 
2000 was about 20 times fainter in the 2-12 keV band than the 
1998 and 1999 events. Assuming similar spectral characteris¬ 
tics for the two events, we conclude that total X-ray fluence 
in the 2-100 keV band was 1.1 erg cm“^ in the interval be¬ 
tween 1998 and 2001. To get the average mass transfer rate, we 
need to know the recurrence time between major outbursts and 
the distance. According to Chen, Shrader, & Livio (1997), the 
sky coverage provided by scanning X-ray instruments has been 
quite good since about 1988, so it seems likely that the recur¬ 
rence time of XTE J1550-564 is at least ten years. The distance 
is not well known, but it is probably about 6 kpc to within a 
factor of two (Table 4). In Table 5 we give total isotropic en¬ 
ergies and average mass transfer rates for assumed distances of 
3, 6 , and 10 kpc and recurrence times of 10 and 50 years (we 
assume an accretion efficiency of 10 % and that the mass trans¬ 
fer rate from the secondary star is the same as the mass transfer 
rate onto the black hole). The mass transfer rates are between 
about 1 X 10 and 7 X 10 yr ', which suggests the star 
is evolving on a relatively long nuclear timescale. 

GRO J1655-40 might be a similar system in terms of its 
evolutionary state. The average mass transfer rate in GRO 
J1655-40 has been estimated to be about 1 x yr“' 

(van Paradijs 1996). However, the secondary star is located in 
the Hertzsprung gap in the Hertzsprung-Russell diagram (e.g. 
Orosz & Bailyn 1997), so one would expect a mass transfer 
rate of about 1 x 10 “^Mq yr“^ (e.g. Kolb et al. 1997). Regos, 
Tout, & Wickramasinghe (1998) have suggested that the sec¬ 
ondary star in GRO J1655-40 is still in the core hydrogen burn¬ 
ing phase, and they have constructed evolutionary models with 
mass transfer rates much closer to the observed value. It re¬ 
mains to be seen if a similar model can be constructed for XTE 
J1550-564. 

Although the presence of a black hole in XTE J1550-564 is 
now firmly established, follow-up observations would be de¬ 
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sirable to confirm and improve upon our results. The statistical 
error in the mass function can be easily reduced by adding more 
radial velocities. Spectroscopic observations with higher reso¬ 
lution should be obtained so that the rotational velocity of the 
companion star can be measured more precisely. Our inclina¬ 
tion limits can be improved by obtaining multicolor light curves 
with much better phase coverage. 


5. SUMMARY 

Our optical spectroscopic observations of the companion star 
in XTE J1550-564 have established the orbital period of Psp = 
1.552 ±0.010 days and the semiamplitude of fG = 349 ± 12 km 
s“'. The optical mass function is f {M) = 6.86 ±0.71 Mq. The 
rotational velocity of the companion star is Vrotsin; = 90± 10 
km s“', which implies a mass ratio of Q= M\/M\ = 

(Icr). This result should be confirmed with better data. Erom 
light curve modelling we find the most likely value of the 
mass of the compact object is 9A1Mq with a Icr range of 
8.36Mq < Ml < 1O.76M0. If we apply our tentative value 
of Vrotsin/ = 90 ± 10 km s“' as an additional constraint in the 
ellipsoidal modelling, we find the most likely value of the 
mass of the compact object is 1O.56M0 with a Icr range of 
9 . 68 M 0 < Ml < 11.58M0. The black hole nature of XTE 
J1550-564 is therefore firmly established by this work. 
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Fig. 1.— (a): The complete YALO V-band light curve of XTE J1550-564 is shown to the left of the arrow located near day 1050. The arrow indicates the time 
of our spectroscopic observations and our F-band photometric observations with VLT/FORSl, which are shown here, (b) The RXTE ASM 2-12 keV light curve of 
XTE J1550-564 (daily averages). They-axis scale is logarithmic to highlight the weaker activity in early 2001 (the measurements with negative flux are not shown). 
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Spectral Type 

Fig. 3. — (a): The rms values of the polynomial fits to the difference spectra as a function of the template spectral type for the templates obtained with VLT/FORSl 
(luminosity class III and IV). The plotting symbols indicate the grism used, (b): Same as (a) but for the templates observed at CTIO. Here the plotting symbols 
indicate the luminosity class of the template. 
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5000 5500 6000 6500 
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Fig. 4.— The results of the spectral decomposition, (a) Top: The dereddened restframe spectrum of XTE J1550-564, smoothed to 5 pixels and offset by 0.5 units, 
(a) Middle: The spectrum of the K3 III star HD 181110, scaled by 0.48 and offset by 0.5 units, (a) Bottom: The difference spectrum, which represents the spectrum 
of the accretion disk, smoothed to 5 pixels with no offset. The interstellar features are denoted by IS and the telluric feature by TEL. (b) Top: The dereddened 
restframe spectrum of XTE J1550-564, smoothed to 5 pixels and offset by 0.2 units, (b) Middle: The spectrum of the G8IV star HD 157931, scaled by 0.74 with no 
offset, (b) Bottom: The difference spectrum smoothed to 5 pixels with no offset. 
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Fig. 5.— (a): The rms values of the polynomial fits to the difference spectra as a function of the input value of the rotational velocity Vi-otsin/. The lowest rms 
value is for Vrotsini = 90 km s“^ for the G8IV template HD 157931 (filled circles) and Viotsinr = 95 km s“^ for the K3in template HD 181110 (filled squares), (b): 
The full width at half maximum of the cross correlation peaks of three different comparison star spectra as a function of the input value of the rotational velocity 
Vrot sink In all cases the template spectrum for the cross correlation analysis was the K3III star HD 181110. The lines show the width of the cross correlation peak 
for the XTE J1550-564 restframe spectrum measured using two different template spectra (the dotted line for the K3III template HD 181110 and the dash-dotted 
line for the G8 IV template HD 157931). The results in (a) and (b) show that the rotational broadening of XTE J1550-564 is on the order of 90 km or less, (c): 
Similar to (a) but for a simulated spectrum constructed with a moderate resolution spectrum 70 km s“^ FWHM), using an input value of Vrot sin i = 90 km s“^. 
(d): Similar to (c) but for the simulated spectrum and the moderate resolution templates. The dotted lines show the full range of the width of the cross correlation 
peak for the simulated object spectrum measured with all six templates, and the dash-dotted lines show the width of the cross correlation peak measured with the 
two best matched templates. The results in (c) and (d) show that the systematic error involved in measuring the rotational velocity in a relatively noisy low resolution 
spectrum is not too large and is perhaps on the order of about 10 km s“^. 
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10 15 5 10 


Fig. 8.— values as a function of the ten fitting parameters for the fitting run A (the run without a Viotsin/ constraint). They are, from left to right and top to 
bottom, the inclination / in degrees, the mass ratio Q, the orbital separation a in solar radii, the orbital period P in days, the phase zero-point To in days since HJD 
2 452 000, the mean temperature of the secondary star Teff, the power-law exponent ^ on the disk temperature profile, the radius of the outer edge of the disk as 
a fraction of the compact object’s Roche lobe radius, the temperature of the inner edge of the disk r^isk in units of 1000 K, and the opening angle of the outer rim 
of the disk in degrees. Fifth or sixth order polynomial fits to the first five curves used to get the confidence regions are shown as the smooth lines. The dashed 
lines denote the 1, 2, and 3cr confidence limits. 
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Fig. 10.— values as a function of various derived parameters for fitting run A (the run without a Vrotsinz constraint). They are the black hole mass Mi in solar 
masses, the secondary star mass M 2 in solar masses, the predicted projected rotational velocity of the secondary star Vrotsinz in km s“^, the radius of the secondary 
star R 2 in solar radii, the surface gravity of the secondary star logg in cgs units, and the bolometric luminosity of the secondary star L 2 in solar units. Fifth or sixth 
order polynomial fits to the curves used to get the confidence regions are shown as the smooth lines. The dashed lines denote the 1, 2, and 3(T confidence limits. 
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Table 1 

Orbital Parameters for XTE J1550-564. 


Parameter 

Value 

Psp, spectroscopic orbital period (days) 

1.552±0.010 

Pph, photometric orbital period (days) 

1.5437 ±0.0005 

K 2 velocity (km s“*) 

349 ±12 

7 velocity (km s“') 

-68 ± 19 

Tq, spectroscopic“(HJD 2,452,000 h-) 

54.296 ±0.014 

To, spectroscopic'’(HJD 2,452,000 h-) 

55.448 ±0.010 

Mass function (Mq) 

6.86±0.71 


Note. — All quoted uncertainties are Icr. 

“The time of the maximum radial velocity of the secondary 
star. 

*’The time of the inferior conjunction of the secondary star. 


Table 2 


Disk Velocities for Black Hole Systems. 


Object Name 

Vd/Kz 

Reference 

H1705-250 

1.11 

1 

GS 1124-683 

1.13 

2 

GRO J0422H-32 

1.17 

3,4 

XTFJ1550-564 

1.20 

5 

A0620-00“ 

1.24 

2 

GRS 1009-45 

1.26 

6 

V404 Cyg 

1.33 

7,8 

GS 2000H-25 

1.35 

9 

A0620-00‘’ 

1.47 

2 


Note. — References: (1) Harlaftis et al. 
1997; (2) Orosz et al. 1994; (3) Orosz & 
Bailyn 1995; (4) Filippenko, Matheson, & 
Ho 1995; (5) This work; (6) Filippenko et 
al. 1999; (7) Casares, Charles, & Naylor 
1992; (8) Gotthelf et al. 1992; (9) Filip¬ 
penko, Matheson, & Barth 1995 

“Observed 1991 profile 

'’Observed 1993 profile 
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Table 3 

Fitted and Derived Astrophysical Parameters for XTE J1550-564 from Light Curve Modelling 


Parameter 

Central 

Value 

(No Vjot sin i constraint) 

Icr range 3 cr range 

Central 

Value 

(Vrotsin/= 90± 10km s ') 

Icr range 3cr range 

Inclination (degrees) 

72.6 

67.0-77.4 

55.2-79.0 

73.5 

70.8-75.4 

58.5-78.4 

Mass ratio 

«21.0 

> 12.0 

>3.6 

8.5 

6.7-11.0 

4.5-22.5 

Orbital separation (/?©) 

11.96 

11.52-12.50 

10.79-15.09 

12.77 

12.35-13.22 

11.53-14.45 

Orbital Period (days) 

1.5435 

1.5430-1.5440 

1.5421-1.5449 

1.5437 

1.5432-1.5441 

1.5421-1.5450 

K 2 (km s-') 

359.4 

350.2-368.6 

331.2-387.4 

361.2 

352.3-370.1 

333.8-388.2 

Black hole mass (M©) 

9.41 

8.36-10.76 

6.81-15.60 

10.56 

9.68-11.58 

8.10-15.15 

Secondary star mass (M©) 

wO.4 

<0.79 

<4.10 

1.31 

0.94-1.64 

0.34-2.73 

Total mass (M©) 

9.62 

8.71-10.86 

7.29-14.83 

11.70 

10.70-12.89 

8.85-16.84 

Secondary star radius (Rq) 

« 1.88 

1.62-2.24 

1.63-4.06 

2.81 

2.55-3.04 

1.93-3.61 

Secondary star luminosity (L©) 

« 1.4 

<2.5 

<6.3 

^4.0 

1.7-5.2 

1.0-7.4 

Secondary star gravity (logg) 

3.46 

3.40-3.56 

3.36-3.84 

3.63 

3.59-3.68 

3.46-3.76 


Table 4 


Distance as a Function oe Mass 


Assumed Secondary 
Mass (M©) 

Radius 

{Rq) 

Mv for 

Teff = 4600 K 

Distance (kpc) 
for Teff = 4600 K 

Icr Distance Range (kpc) 
for 4100 <7;ff< 5100 K 

0.15 

1.47 

5.25 

3.0 

1.4-3.9 

0.50 

2.20 

4.44 

4.4 

2.1-5.7 

1.00 

2.77 

3.97 

5.4 

2.6-7.0 

1.50 

3.17 

3.70 

6.1 

2.9-1.9 

2.00 

3.49 

3.50 

6.7 

3.1-8.7 

3.00 

3.99 

3.23 

7.6 

3.5-9.8 


Table 5 

Total Energy and Mass Transfer Rate 


Assumed Distance 
(kpc) 

Assumed Recurrence 
Time (years) 

2-100 keV Fluence 
(erg) 

Average Mass Transfer 
Rate (M© yr“') 

3 

10 

1.1 X 10‘*5 

6.0 X 10“'° 

3 

50 

1.1 X 10‘*5 

1.2 X 10“'° 

6 

10 

4.3 X 10‘*5 

2.4 X 10“^ 

6 

50 

4.3 X 10‘*5 

4.8 X 10“'° 

10 

10 

1.2 X 10'^'’ 

6.7 X 10“^ 

10 

50 

1.2 X 10'^'’ 

1.3 X 10“^ 




